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20.,  cont. 

The  design  and  initial  development  of  the  stabilized 
laser  gravimeter  system  have  been  previously  described  in 
AFCRL-TR-74-0555.  In  the  design  a mirror  forming  one  end 
of  a Fabry-Perot  optical  cavity  is  mounted  on  a mechanical 
suspension  system.  The  suspension  consists  of  a beam,  end 
loaded  to  obtain  the  required  sensitivity  to  changes  in 
gravitation.  A He-Ne  illuminating  laser  is  locked  t.o  a 
fringe  of  the  Fabry-Perot  cavity.  A He-Ne  laser  stabil- 
ized by  locking  to  a vibration-rotation  absorption  line  of 
a methane  cell  is  used  as  a reference.  The  output  of  the 
system  is  a beat  frequency  between  the  two  lasers  propor- 
tional to  the  gravity  change. 

The  proposed  parametric  excitation  system  applies  a 
periodic  end  load  in  addition  to  the  static  load  on  the 
beam  suspension  system.  It  has  been  shown  theoretically 
and  experimentally  that  significant  (order  of  magnitude) 
gains  in  sensitivity  can  be  realized  by  parametric 
excitation. 
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INTRODUCTION 


The  basic  instrument  used  in  the  tests  discussed  in  this  report 
has  previously  been  described  by  Merchant  [1,2]. 

rhe  instrument  is  a static  gravimeter  as  opposed  to  a pendulum  or 
free-fall  apparatus.  The  gravimeter  consists  of  a mass  elastically 
suspended  so  that  changes  in  gravitation  (or  acceleration'  will  result 
ir.  -i  detectable  relative  motion  of  the  seismic  mass.  The  mass  consists 
of  * mirror  and  holder,  and  the  detection  system  uses  two  lasers.  The 
..lirror  is  the  retro-reflector  of  a high-finesse  Fabry-Pex*ot  cavity,  an 
optical  resonator  with  very  narrow  transmission  fringes.  The  wavelength 
of  an  external  He-Ne  illuminating  laser  is  locked  to  changes  in  the 
cavity  length  by  following  a single  fringe.  The  reference  frequency  is 
provided  by  a second  He-Ne  laser  which  is  stabilized  by  inserting  a 
methane  cell  in  the  cavity,  and  locking  the  laser  frequency  to  the 
molecular  absorption  line  of  methane.  The  output  is  the  beat  frequency 
between  the  stabilized  reference  laser  and  the  cavity-locked  laser.  The 
optical  system  is  an  adaptation  of  the  techniques  used  by  Levine  and 
Hall  [3]  for  a strainmeter. 

The  microgal  resolution  needed  for  a useful  gravimeter  requires 
a very  high  resolution  laser  and  a very  low  natural  frequency  (soft) 
mechanical  suspension  system.  To  obtain  the  necessary  mechanical 
sensitivity,  a suspension  system  consisting  of  an  endloaded  beam  was 
used.  Theoretically,  the  natural  frequency  and  spring  stiffness  of 
the  system  approach  zero  as  the  end  load  approaches  the  buckling  load 
for  the  beam  system  [4,5], 

A vacuum  syster  was  constructed  to  enclose  the  Fabry-Perot  sens- 
ing cavity  and  the  s;  -tem  was  installed  at  a University  of  Washington 
Geophysics  test  site  at  Tumwater,  Washington. 


1.  H.C.  Merchant,  "Stabilized  Laser  Gravimeter,"  AFCRL-TR-74-0355, 

Air  Force  Cambridge  Research  Laboratories,  Hanscom  Air  Force  Base, 
Massachusetts,  July  1974. 

2.  H.C.  Merchant,  E.M.  Hernandez  and  N.D.  McMullen,  "Stabilized  Laser 
Gravimeter,"  Proceedings  of  the  20th  International  Instrumentations 
Symposium,  Albuquerque,  New  Mexico,  May  1974. 

3.  J.  Levine  and  J.L.  Hall,  "Design  and  Operation  of  a Methane  Absorp- 
tion Stabilized  Laser  Strainmeter,"  Journal  of  Geophysical  Research, 
Vol . 77,  No.  14  (May  1972)  pp.  2595-2609. 

4.  L.  Meirovitch,  Analytical  Methods  in  Vibrations,  Macmillan,  New 
York  (1967), 

5.  H.  Lurie,  "Lateral  Vibrations  as  Related  to  Structural  Stability," 
Journal  of  Applied  Physics,  June  1952,  pp.  195-204. 
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The  use  of  parametric  excitation  to  increase  the  sensitivity  of 
the  system  has  been  demonstrated  theoretically  and  on  a laboratory 
model.  The  parametric  excitation  is  introduced  as  a periodic  compo- 
nent superimposed  on  the  static  end  load  of  the  beam  support  system. 

The  sensing  muss  and  mirror  in  the  Fabry- Perot  cavity  are  carried  on 

the  parametrically  excited  support  beam,  resulting  in  an  amplified  mirror 

deflection. 


EXPERIMENTAL  SETUP,  FIELD  TESTS 


The  components  of  the  gravimeter  are  shown  schematically  in  Fig- 
ure 1,  and  their  function  is  identified'.  Figure  2 shows  the  beam  sup- 
port system  and  the  center  mass  containing  a mirror,  which  is  identi- 
fied in  Figure  1 as  the  accelerometer.  Figure  3 shows  the  complete 
Fabry-Perot  cavity  without  the  vacuum  chamber.  The  beam  and  detector 
mass  are  visible  at  the  bottom  in  the  endloading  housing  and  the  He-Ne 
illuminating  laser  for  the  c.  vity  is  visible  at  the  right.  Figure  4 
shows  the  vacuum  housing  for  the  cavity.  The  laser  beam  enters  through 
a quartz  window.  Figure  5 shows  the  Fabry-Perot  cavity  enclosed  in 
its  vacuum  housing,  as  installed  at  the  test  site  at  Tumwater,  Washing- 
ton. The  He-Ne  illuminating  laser  is  visible  at  the  center  and  the 
methane-stabilized  reference  laser  is  in  the  foreground.  The  beat 
frequency  detector  is  visible  at  the  left.  The  details  of  the  con- 
struction and  operation  of  the  system  have  been  reported  previously. 

The  field  tests  were  intended  to  obtain  earth  tides  for  compari- 
son with  predicted  values  over  a time  period  long  enough  to  check  the 
stability  and  sensitivity  of  the  system.  However,  due  primarily  to 
equipment  reliability,  the  maximum  periods  during  which  data  could  be 
obtained  were  hours  rather  than  days  [6] . 

Problems  encountered  in  the  field  tests  were:  degradation  in 

laser  outputs  as  a function  of  time;  multiple  resonances  in  the  Fabry- 
Perot  cavity,  and  multiple-cavity  interaction  between  the  Fabry-Perot 
cavity  and  the  He-Ne  illuminating  laser  cavity;  degradation  of  the 
optical  isolation  system  with  time  (apparently  due  to  moisture  sensi- 
tivity); and  electronic  failures  in  the  servo  System  and  the  long-path 
(Fabry-Perot)  detector  system. 

The  laser  degradation  was  alleviated  by  periodic  laser  realignment. 
The  multiple  resonance  and  cavity  interaction  were  eliminated  by  changes 
in  the  cavity  configuration  and  replacement  of  the  optical  isolation 
system  with  one  with  a higher  transmission  coefficient.  This  higher 
quality  system  also  eliminated  the  initial  problem  with  the  optical 
isolation  s>stem  The  servo  system  was  repaired  before  the  end  of  the 
contract  period,  but  components  to  rebuild  the  long-path  detector  could 
not  be  obtained  in  the  contract  period. 


6.  Laser  Stabilized  Cravimeter  Quarterly  Report,  Nos.  3 through  7, 
Contract  F19628-75-C-0042  (1975-76). 
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Therefore,  as  stated  before,  operation  of  the  system  was  confined 
to  a few  hours.  However,  this  did  allow  the  evaluation  of  the  operating 
parameters  [1]  and  error  sources  as  discussed  in  the  following  section. 


ERROR  ANALYSIS  AND  TESTS 


The  following  sources  of  error  were  examined: 

a)  The  variation  in  vacuum  in  the  Fabry-Perot  cavity. 

b)  The  variation  in  temperature  in  the  Fabry-Perot  cavity. 

c)  Creep  in  the  structural  materials. 

d)  Gravity  noise  due  to  tilt. 

e)  Background  noise  due  to  personnel  movements  and 
local  vehicular  traffic. 

The  vacuum  chamber  surrounding  the  Fabry-Perot  cavity  was  not  construc- 
ted to  obtain  a high  vacuum  but  rather  to  ensure  a stable  condition. 

The  effect  of  a change  in  pressure  is  a change  in  the  speed  of  light 
and  hence  the  cavity  resonance  frequency.  This  frequency  change  is 
reflected  in  an  apparent  gravity  error  (Ag)  as  follows: 


Ag  « (AL)  o)n2 

where 

AT  . 

AL  * - -j-  L 

Af  “ f“*A 


f 


V 

* IT 


f 


A 


» 


and 


(1) 


Ag  * gravity  error 

ojn  = natural  frequency  of  the  sensing  mass  (and  mirror) 

AL  » apparent  cavity  length  change  (deflection  of  sensing  mass) 
L = original  cavity  length 
f = original  cavity  frequency 
Af  = frequency  change  in  cavity 
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f^  = perturbed  cavity  frequency 
V = velocity  of  light  at  the  original  condition 
= velocity  of  light  at  the  perturbed  condition. 


Variation  of  temperature  can  affect  the  syster.  in  two  ways.  The 
first  is  through  a change  in  the  speed  of  light  and  the  second  through 
thermal  deformation  of  the  structure.  The  first  effect  is  handled 
theoretically  in  the  same  manner  as  the  pressure  effect  and  hence  they 
will  be  discussed  together. 

There  have  been  several  investigations,  both  theoretical  and  empiri- 
cal, to  determine  the  relationship  between  the  index  of  refraction  (and 
hence  the  speed  of  light)  and  the  parameters  that  affect  it  [7,  8,  9, 

10].  The  most  recent  is  by  Owens  [10].  His  results  were  used  to  predict 
the  sensitivity  of  the  gravimeter  to  pressure  changes  and  the  aspect  of 
temperature  change  that  affects  the  index  of  refraction. 

The  velocity  of  light  in  a medium  is  related  to  the  index  of 
refraction  by 


where 

C = the  velocity  in  a vacuum 
n * the  index  of  refraction. 

Owens'  results  are: 


7.  H.  Barrell  and  J.E.  Sears,  "The  Refraction  and  Dispersion  of  Air 
for  the  Visible  Spectrum,"  Philosophical  Transactions  of  the  Royal 
Society  of  London,  Series  4,  Vol.  238  (1940)  pp.  1-64. 

8.  Flden  Bengt,  "The  Dispersion  of  Standard  Air,"  Journal  of  the  Optical 
Society  of  America,  Vol.  43,  No.  5 (May  1,  1953),  pp.  339-344. 

9.  Elden  Bengt,  "The  Refractive  Index  of  Air,"  Metrologia,  Vol.  2,  No.  2 
(April  1966)  pp.  71-80. 

10.  James  C.  Owens,  "Optical  Refractive  Index  of  Air:  Dependence  on 

Pressure,  Temperature  and  Composition,"  Applied  Optics,  Vol.  6, 

No.  1 (January  1967)  pp.  51-59. 
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where 


\ • wavelength  (angstroms) 

Ps  ■ partial  pressure  of  air  (millibars) 

Pw  ■ partial  pressure  of  water  vapor  (millibars) 

T ■ temperature  (degrees  Kelvin) 


i_ 

h 

1 

2:>71.34xl0-8 

6487.31  xlO’8 

2 

683939.7x  10"® 

58.058  x 10”8 

3 

130 

-0.7115 x 10"8 

4 

4547.3 x 10’8 

0.08851  x 10"8 

5 

38.9 

3.7  x 10"4 

6 

57.9  x 10"8 

-2.37321 x 10" 

7 

-9,325  xlO”4 

2.23366 

8 

0.25844 

-710.792 

9 

7.75141  x 10’4 

(2) 


The  formula  assumed  atmospheric  air  to  be  made  up  of  dry  carbon  dioxide 
free  air,  water  vapor,  and  carbon  dioxide.  The  partial  pressures  of 
air  aid  water  vapor  required  can  be  determined  by  assuming  ideal  gas 
beha-  ior  at  low  pressure. 

The  assumed  reference  conditions  to  determine  the  effects  of  pressure 
and  temperature  changes  were 
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T = 300°K 
P = 1 ram  Hg 

Relative  Humidity  = 50% 

\ 5 . .SO  1 1 - 33,1)00  A 

The  results  of  applying  the  foregoing  equation  are: 

Pressure 

Ag  * 14.126  milligals/mm  Hg 
Temperature 

Ag  > 16.7  microg&ls/degree  K. 

The  experiment'’!  method  used  to  check  the  pressure  effects  on  the  sys- 
tem was  to  evacuate  the  system,  shut  off  the  vacuum  system,  and  introduce 
a controlled  leak  while  monitoring  the  frequency  of  the  Fabry-Perot 
cavity  as  a function  of  vacuum.  The  test  times  were  short  compared 
to  thermal  delay  times  and  earth  tide  periods.  The  pressure  changes 
selected  were  several  orders  of  magnitude  (2  mm  Hg  or  more)  beyond  the 
system's  operational  fluctuations,  to  obtain  a large  error  output.  The 
predicted  values  for  Ag/mm  Hg  ranged  from  30%  to  40%  below  the  values 
obtained  during  three  pressure  tests.  Therefore,  the  foregoing  calcula- 
tions for  pressure  error  provide  a conservative  estimate  of  the  pressure 
effect. 

To  investigate  the  thermal  effects  on  the  index  of  refraction, 
the  system  was  operated  with  the  vacuum  system  on  and  the  aluminum 
vacuum  housing  heated  at  a series  of  locations.  The  gravimeter  output 
was  monitored  during  the  heating  period  and  for  45  minutes  afterward. 
Temperatures  were  measured  with  thermocouples  attached  to  the  vertical 
quartz  support  rods  (Figure  3,  longest  rods),  to  the  quartz  rods  sup- 
porting the  mechanical  sensing  system  (Figure  3,  short  rods),  and  to 
the  steel  spacers  (Figure  3,  on  top  of  short  quartz  rods),  and  with  two 
thermocouples  suspended  in  the  cavity  next  to  the  laser  beam  path. 

The  other  aspect  of  the  thermal  effect  on  the  gravimeter  is  the 
change  in  cavity  length  due  to  thermal  expansion.  The  support  for  the 
Fabry-Perot  cavity  is  designed  so  that  the  sensor  is  on  re-entry  rods 
of  quartz  and  spacers  of  steel  (Figure  3) . The  support  members  are 
dimensioned  so  that  the  thermal  expansion  of  the  re-entry  supports  cancels 
the  expansion  of  the  lateral  supports  [1,  11).  This,  of  course,  is  the 
case  when  there  are  no  thermal  gradients  between  the  support  elements. 


11.  N.D.  McMullen,  "Methane  Absorption  Stabilized  Laser  Gravimeter: 

Design  of  an  Ultra-Sensitivity  Fabry-Perot  Interferometer  Ac-elerom- 
eter,"  M.S.  Thesis,  Mechanical  Engineering,  University  of  Washington, 
1974. 
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In  the  tests  performed,  thermal  gradients  were  developed  in  the 
supports  that  resulted  in  cavity  length  changes  that  produced  errors 
several  orders  of  magnitude  greater  than  those  that  would  be  predicted 
from  the  temperature  change  (1°K)  obtained  in  the  partial  vacuum.  There- 
fore only  qualitative  statements  can  be  made  about  the  relative  magnitude 
of  the  three  sources  of  error--thermal  expansion,  pressure  changes,  and 
direct  temperature  changes  in  the  cavity.  The  results  indicate  that 
the  direct  temperature  changes  will  cause  errors  that  are  orders  of  magni- 
tude smaller  than  the  pressure  fluctuations.  Therefore  the  most  important 
thermal  concern  is  maintaining  a uniform  temperature  in  the  gravimeter 
structure. 

Creep  is  defined  as  slow  deformation  of  solid  materials  over  extended 
periods  of  time.  Temperature  and  stress  in  the  material  are  parameters 
in  the  creep  process.  The  element  in  the  system  most  sensitive  to  creep 
is  the  support  beam  since  the  material  is  steel  and  it  is  the  most  highly 
stressed.  For  the  creep  to  be  important,  temperatures  on  the  order  of 
one-half  the  melting  temperature  and  stresses  that  are  a significant  per- 
centage of  the  yield  stress  must  occur  [12].  The  maximum  end  load  re- 
quired to  obtain  the  desired  natural  frequency  of  the  sensing  beam  sup- 
port system  (2.5  Hz)  is  14.2  x 106  dynes  with  a resulting  stress  of 
2.8  x 108  dynes/cm2  [13] . The  yield  strength  for  the  beam  material  is 
1.38  x 1010  dynes/cm*.  The  operating  temperature  will  be  less  than  35°C 
while  one-half  the  melting  temperature  is  approximately  750°C.  Therefore, 
no  experimental  creep  tests  were  performed. 

The  effect  of  tilt  at  the  measuring  site  on  the  gravimeter  output 
was  examined  by  introducing  a large  tilt  at  the  base  of  the  instrument 
and  scaling  the  effect  to  geophysical  levels.  The  entire  gravimeter  sys- 
tem was  mounted  on  a granite  slab,  the  edge  of  which  is  visible  in  the 
foregound  of  Figure  5.  A tilt  was  introduced  by  raising  the  edge  of 
the  slab  with  a hydraulic  jack  while  monitoring  the  slab  motion  with 
dial  indicators  supported  on  a frame  over  the  slab.  When  the  results 
were  scaled  to  actual  geophysical  tilt  levels  of  10_9  radian,  the  tilt 
error  was  on  the  order  of  1 ygal  and  therefore  within  the  accuracy  of 
the  gravimeter. 

A discussion  of  the  linear  vibration  analysis  of  the  sensing  system 
is  given  by  McMullen  [11]  and  Bonder  [13].  A discussion  of  the  nonlinear 
aspects  by  Eisinger  [14]  is  summarized  in  the  following  section.  The 


12.  L.  Finnie  and  W.R.  Heller,  Creep  of  Engineering  Materials,  McGraw- 
Hill,  New  York,  New  York  1959. 

13.  G.  Bondor,  "Calibration  and  Error  Analysis  of  a Stabilized  Laser 

Gravimeter,"  M.S.  Thesis,  Mechanical  Engineering,  University  of  Wash- 
ington (in  preparation),  1977.  ' 

14.  K.  Eisinger,  "Parametric  Excitation  of  a Clamped  Beam,"  Ph.D.  Thesis, 
Mechanical  Engineering,  Univeristy  of  Washington  (in  preparation),  1977. 
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purpose  of  the  tests  was  to  identify  operational  problems  empirically. 
Three  types  of  inputs  were  used:  vehicular  traffic  immediately  adjacent 

to  the  test  site  structure  on  the  access  road;  personnel  movement  in  the 
test  site;  and  impacts  adjacent  to  the  test  slab.  The  impact  tests 
were  drop  tests  with  metal-to-metal  impact  to  provide  a higher  frequency 
input  spectrum  than  obtained  by  foot  falls.  An  impact  of  1.9  joules 
(newton  meters)  was  produced  on  the  building  floor  just  adjacent  to  the 
isolation  gap  between  the  pier  and  the  building  floor,  at  2 meters  from 
the  gap,  and  about  3 meters  from  the  gravimeter  on  the  rock  outcropping 
to  which  the  pier  is  attached.  The  results  indicate  that  the  disturbances 
caused  a gravimeter  output  up  to  85  milligals,  which  is  far  above  the 
microgal  sensitivity  required.  However,  the  purpose  of  the  gravimeter 
is  to  measure  phenomena  with  periods  of  hours,  so  such  disturbances 
only  create  spikes  of  noise  in  the  data  which  can  be  filtered  at  the  time 
of  recording  or  during  da  .a  analysis.  The  important  consideration  is 
the  ability  of  the  system  to  maintain  laser  lock  and  provide  continuous 
data.  This  was  accomplished  in  all  the  tests. 


PARAMETRIC  EXCITATION 


The  gravimeter  mechanical  sensor  and  suspension  system  in  Figure  2 
is  shown  schematically  in  Figure  6.  In  the  present  configuration  the 
end  load  [P (t) ] is  a constant,  selected  to  obtain  the  desired  frequency 
and  hence  sensitivity.  By  allowing  the  end  load  to  be  a time-varying 
parameter  with  a specific  frequency  and  mean  load,  it  is  possible  to  obtain 
increased  displacement  of  the  sensing  mass  to  a given  change  in  gravity 
and  hence  an  increased  gravimeter  sensitivity.  This  phenomenon  is  called 
parametric  excitation  or  parametric  amplification  and  has  been  known  for 
many  years.  The  phenomenon  was  observed  by  Faraday  (1831),  Melde  (1859) 
and  Rayleigh  (1883).  Supporting  theory  was  developed  initially  by  Mathieu 
(1868)  and  Hill  (1886).  It  was  later  applied  in  electrical  circuits.  The 
particular  configuration  of  an  endloaded  beam  without  parametric  excita- 
tion has  been  considered  in  a great  many  papers  up  to  the  present  day. 

A detailed  bibliography  on  parametric  excitation  and  related  phenomena 
is  given  in  Appendix  A.  A series  of  papers  by  Rodgers  [15,  16,  17,  18] 
presents  the  development  of  an  electromechanical  parametric  amplifier 
for  use  in  a seismometer. 


15.  P.W.  Rodgers,  "Parametric  Phenomena  as  Applied  to  Vibration  Isolation 
and  Mechanical  Amplifiers,"  Journ.  Sound  and  Vib.,  Vol.  5,  No.  3 
(1967)  p.  486. 

16.  P.W.  Rodgers,  "Sub-Resonant  Response  of  u Mechanical  System,  Para- 
metrically Excited  at  the  Resonant  Frequency,"  Nature,  August  21, 
1965,  p.  853. 

17.  P.W.  Rodgers,  "A  Spring  with  Time-Variable  Stiffness,"  Journ.  of 
Acoust.  Soc.  of  Am. , Vol,  39,  No.  4 (1966)  p.  749. 

18.  P.W.  Rodgers,  "A  Phase  Sensitive  Parametric  Seismometer,"  Bull . of 

Seismological  Soc.  of  Am . , Vol.  56,  No.  4 (1966)  p.  949.  ~ 
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The  analytical  and  experimental  models  were  selected  to  be  compatible 
with  the  existing  gravimeter  configuration.  A schematic  view  of  the  ex- 
perimental setup  is  shown  in  Figure  7,  and  a photograph  in  Figure  8.  A 
beam  of  the  same  dimensions  as  the  gravimeter  sensing  element  with  a center 
mass  (including  the  mass  of  the  accelerometer)  ecfual  to  that  of  the  gravim- 
eter mirror  holder  mass  is  visible  in  the  fixture.  The  static  end  load 
is  provided  by  the  lever  system  visible  at  either  end.  The  dynamic  end 
load  is  provided  by  piezoelectric  elements,  shown  in  Figure  9 as  a cylinder 
at  the  center  of  the  figure  (the  final  configuration  was  a cylindrical 
stack  of  ceramic  disks  for  greater  output) . The  flexure  support  to  main- 
tain the  beam  boundary  conditions  while  allowing  the  end  load  to  be  applied 
is  shown  at  the  left.  The  complete  experimental  setup  is  shown  in  Figure 
10.  The  parametric  system  is  supported  on  shock  cords  and  driven  by  an 
electromagnetic  exciter  through  a direct  current  amplifier. 

The  first  experiments  repeated  Rodgers'  tests  to  check  out  the 
effect  in  the  primary  parametric  amplification  mode  for  the  present 
configuration.  The  results  are  shown  in  Figure  11.  Trace  A is  the  beam 
response  to  an  endload  signal  at  twice  the  beam's  fundamental  frequency. 

The  low  level  response  is  due  to  imperfections  in  the  beam  and  its  align- 
ment in  a gravity  fie.  d.  Trace  B is  the  beam  response  to  a low  level 
base  input  (perpendicular  to  the  plane  of  the  beam)  with  no  variation 
in  the  end  load.  Trace  C is  the  result  of  the  same  base  input  and  periodic 
end  loads  together,  and  demonstrates  parametric  excitation.  In  the  present 
application,  the  system  must  exhibit  parametric  amplification  for  an  input 
signal  with  a frequency  much  less  than  the  natural  frequency  of  the  sensing 
element  which  will  be  identified  here  as  secondary  parametric  excitation. 
The  variables  that  can  be  controlled  are  the  magnitude  of  the  static  end 
load  and  the  frequency  and  magnitude  of  the  dynamic  end  load.  Sinusoidal 
end  loads  were  used  in  the  experimental  and  theoretical  studies.  The 
following  matrix  distinguishes  between  primary  and  secondary  parametric 
excitation. 


Dynamic  endload 

fe) 

Base  motion  v16' 

Phase 

Frequency 

Amplitude 

Frequency 

Amplitude 

Primary 
parametric 
amplification 
(Melde's  effect) 

2 fn 

any 

-fn 

any 

*0° 

Secondary 

parametric 

amplification 

< f 
— n 

F(w2) 

« f 
n 

any 

any 

fn  = system  natural  frequency;  w = end  load  frequency 
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If  the  following  assumptions  are  made  about  the  system  shown  in 
Figure  6, 

a)  Constant  density,  area,  modulus,  and  moment  of  inertia 
along  the  beam 

b)  Small  deflections 

c)  Linear  viscous  damping  for  transverse  motion  (B) 

d)  The  same  input  motion  (acceleration)  at  both  support 
points  (vt|g)) 

e)  End  load  P(t)  = Po  + cos  tot, 

the  resulting  equations  of  motion  are  [14] : 


'tt 


Bvt  + 


kAG  / . \ EA  / \ EA  /.  3\ 

y vx~vxx)  " y (uxvx/x  " 2y  (vx/j 


-v 


(g) 

tt 


(3) 


with  boundary  conditions:  v(Jl,t)  * v(-i,t)  * 0 
, El  , . kAG  / . „ i n 

♦tt  ‘ X^xx  + X"  (*  ■ V * 0 


(4) 


with  boundary  conditions:  ,t)  = 0 


and 


u 


tt 


EA  EA  / 2\ 
y uxx  ” 4y  (vx  jx 


■ 


0 


with  boundary  conditions: 


ux(«,t)  + 


(l,t) 


° 1 

' EA  - EACOSWt 


u(o,t)  « 0, 

where  u,v,  and  <j>  are  the  coordinates  shown  in  Figure  6,  and 

y » beam  mass  per  unit  of  length 
B * damping  coefficient 
A * area 

E * Young's  modulus 
G ® shear  modulus 
I * area  moment  of  inertia 
j * polar  moment  of  inertia 
k = cross  section  shape  factor 
i = beam  length/ 2 
w = end  load  frequency 


(5) 
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A Galerkin  analysis  follows  in  which  a solution  in  the  form  of  a series 
of  comparison  functions  satisfying  the  boundary  conditions  is  assumed. 
For  this  analysis  it  is  convenient  to  combine  the  three  equations.  For 
this  to  be  accomplished,  equation  (5)  must  be  solved  and  the  value  sub- 
stituted in  (3) . Closed  form  solution  of  equation  (5)  is  possible  if 
the  nonlinear  coupling  term  (vx/2)2  is  small  when  compared  to  ux.  This 
depends  upon  the  magnitude  of  the  input,  v-tt^,  and/or  P0  and  Pj.  For 
the  conditions  shown  below  that  are  required  for  parametric  excitation 
in  the  gravimeter,  the  nonlinear  term  may  be  neglected  [14].  The  re- 
sult is: 


EI  vxxxx  - EE  P1  X*4  si"  Xx  cos  Mt  vx 


* 3 EE  P1  X*3  CM  Xx  cos  “l  vxx 


EI 

+ 3 Pj  XX‘  sin  Xx  cos  uit  v 
EI 


. JL  p v 
xxx  kAG  o xxxx 


KAG  P1  xX  cos  Xx  cos  ut  vxxxx 


2 2 2 
, E I 3 - E*I  3 HI  2 

+ 3 TST  vxx  + 9 "KT  vxvxxvxxx  + 7 IT  vx  vxxxx 


- 6 EE  vtxx  ■ “lx)  EE  vttxx  - P1  XX2  sin  Xx  nos  uit  vx 

♦ Po  vxx  * P1  XX  nos  Xx  cos  uit  vxx  • f EA  ^ vxx 

♦ 8vt  . u(x)  vtt  - Jvxxtt  ♦ JJq  p!  “ XX2  sin  Xx  cosuit  vx 

J 2 j 2 

♦ 2 P1  wxX  sinXx  sinwt  vxt  - ^ Pj^  xX  sin  Xx  cos  wt  vxtt 

J 2 J 

- Pj  a)  xX  cos  Xx  cos  wt  v^  - w\X  cos  Xx  sinwt  vxxt 


* BE  Po  vxxtt  * EE  PlXX  cos  Xx  nos«t  »xxtt  - 


, JE  , JE  3 JE  2 q - 

3 EG  vxvxxvxtt  " ^E^vxvxtvxxt  * 2 kGVxVxxtt 


kAG  ttt 


= -Ufx)  V (g)  - U;(x)J  V (g) 

v.  . t- '*■*■+*■  » 


p(x)J 

* kAG"  Vtttt  = vtt'  - TcAG"  vtttt 


(6) 
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subject  to  the  boundary  conditions  (homogeneous  equation): 


v(A,t)  = v(-4,t)  = 0 


vx(4.t)  = <K4.t)  + Y(H,t)  = 0 
vx(-i,t)  = ♦ Y(-M)  = 0 , 


where  the  variables  are  as  defined  previously  and 


- ( 

* \ EA 


2X1/2 


EA7 


and 


The  system's  response  of  interest  is  the  fundamental  (first  symmet- 
ric) mode.  Therefore,  we  use  only  one  term  in  the  approximate  solution, 
with  a comparison  function  approximating  the  first  mode  of  the  linearized 
system  . The  assumed  solution  is: 

v(x,t)  • Cx)q(t)  , 


where  q(t)  describes  the  motion  of  the  beam  center.  This  assumption 
results  in  a fourth-order  equation  in  the  variable  q.  For  a slender 
beam  of  the  gravimeter  configuration  and  for  a low  excitation  (base 
motion)  frequency,  the  equation  can  be  approximated  by  a damped,  forced, 
nonlinear  Mathieu  equation, 


a2iq 


aaq 


♦ (S, 


01 


~ 3 ~ 

♦ aQ2  cos  ujt)  q+  Bflq  + Bj 


2 ~ >2 
qq  + 2Bjq  q 


(7) 


The  equation  was  numerically  integrated  and  a linearized  version  of  the 
equation  was  analyzed  for  conditions  consistent  with  the  gravimeter. 

The  results  were  within  a percent  of  each  other.  The  linearized  equa- 
tion, which  was  used  to  prepare  the  stability  and  amplification  plots 
shown  in  Figures  13-15,  is: 

2 

-*•  2£cx^f-  + (a2  -Yjcoswt)  q ■ -YjA  sinwt  , (8) 

dt 


where 
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a * a0l/52l 


^1  ~ ”®02^®21  * ^ s acc*  amplitude  in  m/sec^ 


a21  Va21 


“21  ’ 3“o  •«f*w(i)1k**-rl 


-J  (if  . J p /»» 


' w - ire  ro  U/ 

“12  * ar  *re  pi  “x*  x4  * 2 eg-  ^ «x  '^(J)2 

^ « (?r p.  lap  * m4i 

^■pi  *\> 

■ ro  x*(f)4  * *r  xw4  - xx(J)2 
• w ire  “2  X2X4  * 5^  o,2  xx(f)2j . 

7Ms  equation  can  be  transformed  to  a standard  form, 
d2y 

— <r  + (a  - 2q  cos  2z)  y • F(z) 
dz 


where 


F(z)  « B sinfiz 


I 


q(z)  » e‘<Z  y(z)  , u>t  * 2z 

as4\(l-(2)»a-K2 


q * 2 


Y,A 

B = -4  ~ eKZ 

<2/ 
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a = 2 - 

U) 


K = 2^  • 


fhe  solution  for  the  homogeneous  equation  is  given  in  terms  of 
Mathieu  functions  (integral  and  fractional  order) . A stability  plot  for 
the  equation  as  a function  of  the  parameters  a and  q is  shown  in  Figure  12. 

The  solutions  of  the  forced  equation  show  that  the  desired  response  at 
the  frequency  of  the  input  is  modulated  by  other,  multiple  (higher)  fre- 
quency components.  An  analysis  of  the  solution  in  the  stable  regions 
of  Figure  12  shows  that  both  parametric  amplification  and  parametric 
attenuation  can  occur.  Figures  13  and  14  show  the  amplification  (or 
attenuation)  curves  for  the  first  three  regions.  Figure  15  shows  the 
effect  of  damping  in  regions  2 and  3.  Figure  16  gives  the  computer  cal- 
culated response  for  parameters  in  region  3 (amplification)  showing  the 
modulated  signal.  Figure  17  shows  the  low-pass  filtered  signal. 

The  above  data  can  be  used  to  select  the  coefficients  of  the  Mathieu 
equation  that  will  result  in  amplified  responses  tha . correspond  to  the 
stable  region  in  Figure  12.  A more  useful  presentation  for  design  purposes, 
however,  is  given  in  Figures  18  and  19. 

If  a constant  value  for  the  parameter  q in  Figure  12  is  considered, 
it  is  possible  to  plot  the  points  of  intersection  with  the  characteristic 
curves  (aj  and  bj)  in  Figure  12  in  terms  of  the  static  end  load  (P0) , 
the  total  dynamic  endload  magnitude  (P0  * Pi) , the  endload  frequency 
(w) , and  the  transverse  frequency  of  the  beam  (a) . One  example  of  such 
a plot  is  shown  for  q ■ 1 in  Figure  18,  which  gives  the  static  end  load 
(normalized  by  the  critical  buckling  load,  Pcr)  as  a function  of  the  end- 
load frequency  (normalized  by  the  transverse  undamped  natural  frequency  of 
the  r.onendloaded  beam,  ap)  Figure  19  is  another  example  for  q * 1 that 
plots  the  magnitude  of  '„ne  dynamic  end  load  as  a function  of  the  trans- 
verse beam  frequency.  A set  of  curves  (with  a range  of  q values)  would 
allow  the  selection  of  the  static  end  load,  dynamic  end  load,  and  endload 
frequency  to  attain  a desired  point  in  the  stability  and  amplification 
regions  shown  in  Figures  12  through  14. 

The  results  of  one  of  the  experiments  conducted  to  verify  the  fore- 
going theory  of  secondary  parametric  amplification  are  shown  in  Figure  20. 

The  system  parameters  are  shown  on  the  figure;  the  experimental  setup 

was  the  same  as  shown  in  Figures  7-10.  The  high  frequency  modulation  is 

also  apparent  in  the  case  of  secondary  amplification.  ' ’■ 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  approach  taken  to  develop  a gravimeter  using  a mechanical 
sensor  in  the  form  of  an  endloaded  beam  and  a stabilized  laser  sensing 
system  still  holds  promise.  No  experimental  or  theoretical  reasons  have 
been  identified  that  would  prevent  the  operation  of  such  a gravimeter. 
However,  it  has  not  been  possible  within  the  contract  period  to  operate 
the  instrument  at  its  theoretical  sensitivity  for  a period  of  time  long 
enough  to  obtain  earth  tides.  The  potential  of  using  parametric  excita- 
tion to  further  increase  the  sensitivity  of  this  gravimeter  or  for  a simi- 
lar application  has  also  been  demonstrated  theoretically  and  experimentally. 

0 

It  is  recommended  that  a system  using  visible  red  light  (6328  A) 
and  an  iodine  rather  than  a methane  absorption  line  be  constructed.  This 
would  produce  the  following  benefits: 

a)  A decrease  in  cavity  size,  with  a corresponding  decrease  in 
potential  sources  of  errors 

b)  A direct  current  rather  than  RF-excited  system,  which  would 
alleviate  noise  problems  associated  with  the  detection  system 
in  the  present  configuration 

c)  The  direct  alignment  of  the  system  with  the  operating  red 
light  beam  rather  than  with  an  auxiliary  laser  system, 
which  would  allow  more  accurate  initial  alignment  and 
routine  monitoring  of  the  system. 

It  is  not  recommended  that  parametric  excitation  of  a gravimeter 
system  be  carried  beyond  its  present  state  of  development  until  either 
the  present  methane-stabilized  system  or  an  iodine-stabilized  system  has 
successfully  demonstrated  its  ability  to  detect  earth  tides. 
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Figure  3.  Fabry-Perot  sensing  cavity. 
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Figure  4.  Cavity  vacuum  chamber 
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Figure  6.  Schematic  of  mechanical  sensor. 


— OSCILLOSCOPES  — 


1.  CLAMPED  BEAM  4. 
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PIEZOELECTRIC  CRYSTALS 
((7)  PZT-4) 

STATIC  PRELOAD  LEVEL 

ELECTROMAGNETIC  EXCITER 
(LING  420) 


Figure  7.  Schematic  of  test  setup. 
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Figure  9.  Excitation  system  components. 
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Figure  10.  Parametric  excitation  system. 
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BEAM  RESPONSE 


BASE  INPUT  MOTION 
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NATURAL  FREQUENCY;  NOMINAL  BASE 
MOTION  INPUT;  PHASE  LOCKED 


Figure  11.  Primary  parametric  amplification  phenomenon, 
clamped-c lamped. 
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CZZ]  UNSTABLE 
I 1 STABLE 

3 = CONSTANT  LINES 

fi  = CONSTANT  LINES 

(P>  NUMBERING  SYSTEM  FOR  STABLE  REGIONS 


l;ij»iire  12.  Mathieu  equation  stability. 
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Figure  14.  Secondary  parametric  amplification 
characteristic  region  3,  Figure  12 


Figure  If).  Secondary  parametric  amplification  with  damping. 
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MID-BEAM  DISPLACEMENT,  X 10  cm 


TIME,  SECONDS 


Figure  16.  Numerical  analysis  of  a clamped  straight  beam. 


(D>  DISTORTION  DUE  TO  DATA  REDUTION  PROCESS 


Figure  17.  Low-pass  filtered  response  of  a clamped-clamped  beam. 


Figure  18.  Interaction  of  system  parameters  a,  q,  PQ  and  w 
for  q = 1.0,  Figure  12. 
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ENDLOAD  MAGNITUDE,  P0i  P| 


UNDAMPED  NATURAL  FREQUENCY  IN 
TRANSVERSE  BEAM  AXIS,  a/a0 


Figure  19.  Interaction  of  system  parameters  a,  q,  P , P and  a 
for  q = 1.0,  Figure  12.  0 
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0123456789  10 
TIME,  SECONDS 

fa)  NO  DYNAMIC  END  LOAD 


0.096  V 


0123456789  10 
TIME,  SECONDS 


(b)  DYNAMIC  END  LOAD 


SYSTEM'S  NATURAL  FREQUENCY  ■ 9 Hz 
ENDLOAD  FREQUENCY  - 8.7  Hz 

AMPLITUDE  ■ 2.2 x 10s  DYNES 
BASE  FREQUENCY  » 0.33  Hz 


Figure  20.  Secondary  parametric  amplification  phenomenon  for  a 
claaped-clanped  beam.  * 
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APPENDIX  A 


BACKGROUND  AND  LITERATURE  SURVEY,  PARAMETRIC  EXCITATION 


The  phenomenon  of  parametric  excitation  of  oscillatory  systems  has 
been  known  for  many  years.  The  literature  describing  and  evaluating 
this  phenomenon  goes  back  to  the  nineteenth  century.  In  1831  M.  Faraday 

[1]  observed  parametrically  induced  oscillations  on  vibrating  surfaces. 
Later,  similar  discoveries  were  made  and  reported  by  F.  Melde  in  1859 

[2]  and  by  Lord  Rayleigh  in  1883  [3j.  The  mathematics  to  allow  analyti- 
cal treatment  of  the  classical  problems  was  developed  simultaneously  by 
Mathieu  in  1868  [4]  and  Hill  in  1886  [5] . Thus  the  phenomenon  of  para- 
metric excitation  was  first  discovered  on  mechanical  systems.  Years 
later,  however,  these  effects  were  studied  and  applied  effectively  to 
electrical  circuits  by  Brillouin  in  1897  [6]  and  by  Poincare  in  1907  [7], 
In  the  years  to  follow,  the  phenomenon  of  parametric  excitation  of  os- 
cillatory systems  has  found  an  ever  increasing  number  of  applications, 
many  of  which  are  in  the  field  of  electronic  circuit  theory. 

Investigations  of  transverse  vibrations  of  uniform  beams  have  been 
conducted  since  the  beginning  of  this  century.  Early  investigators 
were  mainly  concerned  with  the  general  form  of  the  governing  differential 
equation  of  motion  for  free  transverse  vibration  of  uniform  beams  for 
various  end  conditions.  Second-order  corrections,  such  as  terms  for 
transverse  shear,  rotary  inertia,  damping,  etc.,  were  also  proposed. 

Most  of  this  early  work  can  be  found  in  papers  by  Holzer  [8],  Timoshenko 
[9],  Muto  [10],  Sezawa  [11],  and  Suyehiro  [12].  The  best  known  result  of 
these  studies  is  the  familiar  Timoshenko  Beam  Equation,  which  includes 
both  transverse  shear  and  rotary  inertia  of  the  differential  beam  sec- 
tion (sec  also  [31]).  The  effects  of  compressive  end  loads  on  a slender 
column,  in  particular  on  its  frequency  behavior,  were  initially  studied 
by  Sezawa  [13,  14],  Howland  [15],  Cowley  [16],  and  Levy  [17], 


More  recent  investigators  of  transverse  vibrations  of  uniform  beams 
arc  again  concerned  with  second-order  nonlinearities  in  the  beam  equation. 
In  general,  such  nonlinearities  arise  from: 


•%  \ 


b) 

c) 

d) 

e) 


longitudinal  constraints  which  arc  either 
immovable 

some  initial  curvature  of  the  beam 
inclusion  of  transverse  shear 
rotational  inertia  terms 
longitudinal  inertia  terms. 


movable 


or 


Some  analyses  that  focus  primarily  on  the  nonlinear  elasticity  terms, 
such  as  a)  to  c) , were  published  by  Eringen  [18],  Woinowski-Krieger  [19], 
McDonald  and  Raleigh  [2],  Burgreen  [21],  Wahi  [22],  Eisley  [23],  Bennett 
and  Eisley  [24] , Tseng  and  Dugundji  [25,  26],  Srinivasan  [27]  and  Lurie 
[28J . The  resulting  equations  of  motion  obtained  by  the  above  authors 


are  restricted  to  vibration  amplitudes  that  are  less  than  or  equal  to 
the  thickness  of  the  beam  as  well  as  to  fixed  end  coordinates  (as  com- 
pared to  constant  end  loads) . Solutions  to  these  equations  have  been 
obtained  as  follows:  by  Eringen  [18]  by  means  of  perturbation  techniques; 

by  McDonald  and  Raleigh  [20]  using  an  expansion  into  Jacobian  elliptic 
functions  of  time;  by  Burgreen  [21]  by  assuming  a starting  mode  shape 
and  solution  of  the  resulting  integral  equation  of  motion  using  elliptic 
integral  representation  of  time;  by  Bennett  and  Eisley  [24]  using  the 
forced  response  in  a three-mode  case;  by  Tseng  and  Dugundji  [25,  26] 
using  Galerkin’s  method  and  a solution  of  the  resulting  Duffing  equation 
by  the  harmonic  balance  method;  and  by  Srinivasan  [27]  using  the  Ritz 
averaging  method  in  both  space  and  time  variables  to  study  the  steady 
state  forced  response  in  two  modes.  The  papers  of  Tseng  and  Dugundji 
are  of  particular  interest  since  in  addition  to  simple  harmonic  motion 
they  address  sub-  and  superharmonic  motion.  Other  papers  that  also 
address  this  subject  were  published  by  Snowdon  [34,  35]  and  Sridhar, 

Nayfeh  and  Mook  [36] . Using  the  Timoshenko  Beam  Equation,  which  con- 
siders transverse  shear  and  rotary  inertia  terms  non-negligible,  Snowdon 
computes  displacement  amplification  curves  for  various  end  conditions 
of  free  vibration  in  his  first  paper  [34]  and  of  forced  vibration  in  his 
second  paper  [35],  The  paper  by  Sridhar,  Nayfeh  and  Mook  [36]  departs 
from  a general  second  order  nonlinear  equation  into  which  the  nonlinear 
beam  equation  always  can  be  cast  and  presents  solutions  using  perturba- 
tion techniques,  i.e.,  the  method  of  multiple  scales  as  developed  by 
Nayfeh  (see  [58])  and  other  authors.  Subharmonics,  superharmonics  and 
combination  resonances  are  investigated;  the  study  considers  a multiple 
degree  of  freedom  system. 

In  contrast  to  the  papers  mentioned  previously,  in  which  primarily 
"elastic  nonlinearities"  are  considered,  Atluri  [29]  addresses  the  prob- 
lem of  modification  of  the  equation  of  motion  by  inclusion  of  longitudi- 
nal and  rotational  inertia  terms  only.  Using  Galerkin's  method  to  derive 
the  equations  of  motion,  he  finds  solutions  to  these  equations  by  means 
of  perturbation  techniques.  1 

Another  kind  of  continuous  beam  problem  is  the  one  of  constant  end 
loads  as  compared  to  constant  end  coordinates  which. was  considered  above. 
Such  problems  have  been  discussed  extensively  by  Lord  Rayleigh  [30]  and 
Timoshenko  [31].  While  all  previous  references  are  restricted  to  ampli- 
tudes of  oscillation  that  do  not  exceed  the  thickness  of  the  beam,  a 
paper  by  Woodall  [32]  addresses  large  amplitude  oscillations  of  thin 
elastic  beams.  The  solutions  to  his  equations  of  motion--which  he  derives 
by  means  of  the  Galerkin  method--use  the  perturbation  method  and  the  finite 
difference  method.  The  influence  of  a concentrated  mass  on  the  free  vibra- 
tion of  a uniform  beam  is  discussed  in  a paper  by  Maltbaek  [33],  Starting 
with  a sixth-order  frequency  determinant,  his  paper  lists  expressions  for 
the  beam  frequencies  according  to  end  conditions  and  location  of  the 
concentrated  mass. 

Another  class  of  literature  is  the  one  that  describes  the  application 
of  the  phenomenon  of  parametric  excitation  to  certain  mechanical  systems. 

As  noted  previously,  parametric  amplifiers  have  (although  they  were  dis- 
covered in  mechanical  systems)  found  in  recent  years  increasingly  more 

\ 
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applications  in  electrical  engineering.  Some  of  the  more  recent  classi- 
cal references  in  this  field  arc  by  Blackwell  and  Kotzebue  (37]  and  Tucker 
[38].  Application  of  the  parametric  phenomenon  to  mechanical  oscillatory 
systems  often  leads  to  the  classical  problem  of  the  dynamic  stability 
of  slender  columns  that  are  subjected  to  periodic  longitudinal  forces.  The 
basic  concern  that  led  to  most  studies  in  this  area  is  the  occurrence  of 
excitation  of  the  fundamental  transverse  beam  frequency  at  longitudinal 
forcing  frequencies  other  than  the  fundamental  frequency.  This  presents 
a serious  design  problem,  in  particular,  for  some  aircraft  structures.  Two 
of  the  early  authors  on  this  subject  were  Utida  and  Sezawa  [39].  In  their 
report,  the  authors  discuss  extensive  testing  of  a brass  cantilever  beam 
with  a concentrated  mass  at  its  end;  the  periodic  end  load  was  applied 
through  a magnetic  field  at  the  free  end.  The  mathematical  treatment  of 
the  problem  was  limited  to  stability  investigations  of  the  well  known 
Mathieu  equation  into  which  the  problem  was  cast.  It  is  interesting  to 
note  that  the  authors  found  that  instabilities  of  the  system  (or  conversely 
parametric  excitation)  occur  at  ratios  of  the  natural  transverse  beam 
frequency  of  1/2,  2/2,  3/2,  4/2,  5/2,...,  with  a decrease  of  the  resonant 
amplitude  with  increase  of  this  ratio.  Similar  results  have  been  found  by 
many  authors  to  follow. 

A somewhat  different  analytical  approach  to  this  problem  was  taken 
by  Mettler  [40,  41,  42],  In  his  works,  Mettler  used  energy  principles, 
in  particular  the  principle  of  virtual  displacement  (see  also  Marguerre 
[43])  and  Hamilton's  Principle  to  derive  the  equations  of  motion  of  a 
pinned,  uniform,  homogeneous  beam  under  periodic  axial  loading.  In  his 
papers,  Mettler  reduces  the  equations  of  motion  to  a system  of  n homogen- 
eous, linear,  second-order  differential  equations  with  periodic  coeffi- 
cients in  the  form  of  a general  Mathieu  equation' which  he  solves  with 
a so-called  double  series  development.  One  is  a trigonometric  series 
and  the  other  a perturbation  series.  He  obtains  plots  of  instability 
regions  that  are  similar  to  those  obtained  by  Utida  and  Sezawa  [39]  but 
they  are  developed  by  purely  analytical  methods.  These  papers  are  par- 
ticularly noteworthy  since  they  depart  from  the  basic  principles  of 
the  theory  of  elasticity  to  derive  the  required  pnergy  relations.  Investi- 
gations based  on  these  fundamental  papers  were  published  by  Mettler  and 
Weidenhammer  [44]  and  Wcidenhammer  [45,  46,  47].  An  interesting  finding 
reported  in  the  paper  by  Mettler  and  Weidenhammer  [44]  is,  that  if  a beam 
is  subjected  to  harmonic  end  coordinate  displacements  the  spring  behavior 
is  of  the  hardening  type;  if  the  beam  is  subjected* to  harmonic  end  loads 
only,  it  essentially  behaves  linearly  with  theoretically  unlimited 
growth  of  amplitude  (neglecting  second-order  nonlinearities  in  the  equa- 
tion) . If  a concentrated  mass  is  located  at  one  end  and  the  beam  is 
subjected  to  harmonic  end  loads,  the  spring  behavior  is  of  the  softening 
type.  These  analyses  address  a pinned  configuration.  A paper  by  Weidenhammer 
[45]  studies  the  problem  of  a clamped  beam  configuration  with  a periodic 
axial  force.  The  treatment  of  this  problem  is  much  like  that  of  the 
pinned  problems  by  Mettler  [40,  41,  42]  and  in  fact  is  based  upon  the 
same  theoretical  aspects. 

Another  analysis  approach  to  problems  on  flutter  and  autoparametric 
resonances  is  taken  by  Herrmann  and  Hauger  [48]  by  presenting  solutions 
in  the  form  of  Fourier  series.  The  special  problem  of  "snap-through" 


(symmetrical)  and  one-sided  (unsymmetrical)  vibrations  of  columns  is  pre- 
sented in  a paper  by  Min  and  Eisley  [49],  Based  on  a paper  by  Wahi  [22] 
that  assumes  the  space  and  time  variables  are  separable,  the  authors  use 
a normal  mode  solution  in  the  form  of  an  assumed  series.  Limiting  the 
discussion  to  the  assumption  of  two  modes,  the  authors  find  an  approxi- 
mate solution  to  the  nonlinear  equation  of  motion  using  the  harmonic 
balance  method.  Analog  computer  solutions  are  also  used  for  comparison 
with  the  analytical  results. 

Finally,  a series  of  papers  was  written  by  Rodgers  [50,  51,  52,  53]  on 
the  development  of  a mechanical  parametric  amplifier  and  its  application 
as  a phase  sensitive  parametric  seismometer.  In  contrast  to  the  pre- 
vious references,  tuese  papers  discuss  the  general  behavior  that  is  to 
be  expected  of  a damped  single  degree  of  freedom  oscillator  with  a spring 
of  time-variable  stiffness  [52] . The  spring  used  by  Rodgers  [52]  is 
of  the  magneto-mechanical  type.  A noteworthy  aspect  of  this  paper  is  a 
clear  demonstration  of  the  phase  sensitivity  of  the  ratio  of  the  input 
to  the  parametric  response,  e.g.,  the  amplification  curve.  Only  primary 
effects  at  endload  frequencies  of  twice  the  oscillator  natural  frequency 
and  of  the  natural  frequency  and  no  nonlinearities  are  considered. 

Besides  the  above  papers,  a series  of  reference  texts  have  appeared 
that  address  chapters  to  the  phenomena  of  parametric  amplification  of 
oscillatory  motion.  The  references  of  interest,  which  are  by  Minorsky 
[54] , Bolotin  [55] , and  Kononenko  [56] , were  all  originally  Russian 
publications. 

Information  on  the  type  of  data  that  is  to  be  measured  by  a gravim- 
eter, which  is  the  proposed  application  of  this  device,  is  given  in  a 
text  on  earth  tides  by  Melchior  [57], 
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